Effects of diffusion and particle size in a kinetic model of catalyzed reactions 
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We study a model for unimolecular reaction in a supported catalyst including reactant diffusion 
and desorption, using analytical methods and scaling concepts. For rapid reactions, enhancing 
surface diffusion or increasing particle size favors the flux of reactants to the catalyst particles, 
which increases the turnover frequency (TOF). The reactant flux towards the support becomes 
dominant when the ratio of diffusion lengths in the catalyst and in the support exceeds a critical 
value. A peak in the TOF is obtained for temperature-dependent rates if desorption energy in the 
support (Ed) exceeds those of diffusion (Ed) and reaction (Er). Signiflcant dependence on particle 
size is observed when the gaps between those energies are small, with small particles giving higher 
TOF. Slow reactions {Er > Ed) give TOF monotonically increasing with temperature, with higher 
reactant losses in small particles. The scaling concepts can be extended to interpret experimental 
data and results of more complex models. 
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INTRODUCTION 



Modeling heterogeneous catalytic processes is an essen- 
tial tool for catalyst design and for the improvement of 
operating conditions [Ij, i^, |3i, tJ] . Hierarchical approaches 
have to be adopted due to the need of information on a 
wide range of length and time scales, from the electronic 
structure to the reactor design [i| . An important step of 
this approach is the microkinetic modeling, where micro- 
scopic processes such as reaction, diffusion, aggregation, 
and desorption are described by stochastic rules, provid- 
ing information on the efficiency of the catalytic process 
in length scales ranging from a few nanometers to several 
micrometers. 

An important problem to be addressed with these 
methods is the effect of diffusion of reactants through the 
interface between the catalyst particle and the support. 
Several recent experimental papers illustrate the se p he- 
nomena in catalyzed reactions [a, S 0, B B [Hlj 112 
as well as in related problems, such as gas adsorption 
[13, [Hi [Si I where the same materials may be used. 
Morphological features of the catalyst and support and 
physicochemical conditions of operation determine the 
role of diffusion on the performance of the process [3, [13] ■ 
Sometimes these phenomena are called capture-zone ef- 
fects, since a certain region of the support surrounding 
the particle may increase the effective area for capturing 
reactants from the gas phase, but there are also cases 
where a net flux to the support is observed. When only 
one phase can adsorb the reactant from the gas, the terms 
spillover (reactant flux from the catalyst to the support) 
or back spillover (the opposite flow) are used [3 > 
though recently many authors have extended these terms 
to reactants adsorbing in both phases [^ (and here we 
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will use them to facilitate the discussion of the results). 

Due to the large interest in industry, some mod- 
els which incorporated the effects of reactant diffusion 
through the catalyst-support interface were designed for 
certain applications. The simplest models are based on 
rate equations (mean-fleld models) that do not account 
for the spatial heterogeneity of the media where reac- 
tions take place [13, [HI [H, , or that use some type of 
approximation to represent that heterogeneity [24 . [25| . 
Other models represent it through distributions of cat- 
alytic sites in lattices. Most of them are designed to de- 
scribe CO oxidation in different catalysts and conditions 
[H, [13,111, mi , and simplify diffusion and adsorption pro- 
cesses of some species (although other applications have 
also been proposed [la, [13, [IS, [3l|). In order to get a 
deeper insight on the effects of reactant diffusion, Cwiklik 
et al [s?, ^33\ recently simulated simple reaction-diffusion 
models in surfaces with catalytic stripes and squares, as 
well as random distributions of catalytic sites. In cer- 
tain ranges of parameters, they showed monotonic de- 
pendences of the turnover frequency (TOF) on diffusion 
coefficients and reaction rates [s^ . 

A number of other papers aim at a full investigation of 
simple models of reaction and diffusion. In the present 
context, relevant examples are models in lattices with dis- 
tributions of catalytic sites, i. e. with some type of non- 
homogeneity of catalytic activity [H, [H, [H, [13 [H, [1^ . 
Even adsorption-desorption models without surface dif- 
fusion show that the correlations in catalyst particle po- 
sition have nontrivial effects on the TOF, independently 
of adsorbate interactions [s^, [10, [4l[ . Moreover, models 
including diffusion in heterogeneous surfaces show that 
the structures that maximize the efficiency of a catalytic 
process are highly dependent on the rates of the main 
microscopic processes [39| . 

In the present paper, we will propose a one-dimensional 
lattice model for unimolecular reactions in a supported 
catalyst, with reactant diffusion and desorption both in 
the support and in the catalyst particles. Our aim is 
to understand the interplay between these physicochemi- 
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cal mechanisms, particle size and catalyst coverage. The 
model geometry is equivalent to that of Ref. [S^], but 
here we will obtain an analytic solution that facilitates 
the illustration of different possible outcomes. We will 
use scaling concepts to explain the model results, so that 
this framework can be extended to more complex models 
and applications to real systems, where numerical solu- 
tions are usually necessary. Indeed, scaling approaches 
were already shown to be very useful to understand qual- 
itative trends in reaction-diffusion models [i^ . 

Among our results, we will distinguish conditions to 
enhance the net reactant flux from the catalyst to the 
support or vice-versa by varying one of the microscopic 
rates, and we will discuss the effect of increasing the cat- 
alyst particle size. The increase or decrease of the TOF 
will be shown to depend on the relation between diffusion 
lengths of reactants in the catalyst and in the support. 
Some of these results reinforce findings of previous works 
[27I [28I . [3^ . Moreover, under reasonable assumptions 
for temperature-dependent rates, we will show that a re- 
markable increase in the TOF can be obtained if reverse 
spillover regularly fills the catalytic particles with the re- 
actants adsorbed in the support. This feature may be ob- 
served in a large temperature range, a possibility which 
is interesting for applications. The identification of these 
scenarios is possible because the model accounts for the 
inhomogeneity of the catalytic system and, consequently, 
predicts inhomogeneous distributions of reactants, which 
advances over the mean-field models. 



II. THE MODEL 



The physicochemical processes involved in the model 
are illustrated in Fig. 1. 
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FIG. 1; Schematic representation of the model for the supported 
catalyst, with the rates of the physicochemical processes. Non- 
adsorbed reactant A (equation |3]l is represented by white circles, 
while adsorbed reactant A is represented by black circles. The x 
axis used in the analytic solution are shown (for the particle and 
for the support). 

The catalyst particles are represented by segments of 
I sites in a line, separated by d support sites. Assuming 
that a is the size of the lattice site, this corresponds to 
particles of diameter la separated by a distance da. The 
fraction of the support covered by the catalytic material 
is 



This lattice structure may be a reasonable description of 
some model catalysts O , as well as a good approxi- 
mation to the morphology of catalytic clusters supported 
in materials with long and narrow pores [4^1 . 

The particle size is determined by the physico-chemical 
conditions in which the catalytic material is deposited on 
the support (45j . The fraction of the support covered, e, 
is related to the amount of material used to produce a 
sample as well as the particle shape, which determines 
the surface to volume ratio. In oxide supported Pt or 
Pd catalysts, the particle diameter usually vary from 1 
to 50 nm - depending on the growth conditions and if 
sintering occurs, sizes of 100 nm or more can be found. 
Since a is of order of a few angstroms, this typically cor- 
responds to I ranging from 3 to 150 lattice sites. The 
spacing between the particles has a broad distribution 
in catalysts supported in porous materials, but in model 
catalysts they are nearly uniform, usually in the range 
50 — 200 nm [d typically between 100 and 500 lattice 
spacings for most oxide supports). 

The flux of a single chemical species (reactant) towards 
the surface occurs with rate _F, which is defined as the 
number of incoming molecules per site per unit time. In 
most of this work, F will be used to define the time scale 
of the model, so that other frequencies are calculated 
relatively to this quantity. The reactant adsorbs in the 
(randomly chosen) site of incidence if it is empty, either in 
the catalyst or in the support, otherwise the adsorption 
attempt is rejected. Sticking coefficients are set equal to 
1 , since the effects of different values in the catalyst and 
in the support can be incorporated in the corresponding 
desorption rates: and k% in the support and in the 
catalyst particles, respectively (each rate corresponds to 
number of events per site per unit time). Interaction be- 
tween the adsorbates is restricted to the excluded volume 
condition. 

Adsorbed reactants diffuse with coefficient D, which 
for simplicity is assumed to be the same in the support 
and in the catalyst. It means that each reactant attempts 
to execute ID ja^ random steps to nearest neighbor sites 
per unit time. We are assuming that the activation en- 
ergy for diffusion is the same in those regions, which is 
certainly not true for a real catalyst. However, it is not 
a very restrictive assumption for our study because the 
results are interpreted in terms of diffusion lengths and 
scaling ideas are emphasized, allowing a direct extension 
to cases of different D in different regions. 

Since reactants are always in contact with the support 
or the catalyst, this is a surface diffusion model, which is 
a reasonable assumption on metal particles. On the other 
hand, Knudsen diffusion may be more realistic for the 
movement of reactants inside a catalyst pore. However, 
this does not represent a limitation because a suitable 
value of D can be chosen, and interpretations based on 
diffusion lengths are still valid. 

The adopted reaction mechanism for reactant A and 
product B is 

A-ads ^ Bgas, (2) 
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with reaction rate kr assumed to be uniform in the cat- 
alytic region (this rate corresponds to number of events 
per site per unit time). The unimolecular reaction ([2|) 
may represent an Eley-Rideal (ER) mechanism, in which 
the adsorbed species reacts with another species C arriv- 
ing from the gas phase and forms a volatile product: 

-^ads ~^ Cgas ^ Bgas- (3) 

For instance, an application to CO oxidation is discussed 
in Refs. [IHH^I- In this case, kr not only accounts for 
activation of the adsorbed reactant but also for the flux 
of the other reactant from the gas phase. Certainly the 
assumption that k^. is constant in the whole catalytic re- 
gion is not realistic because it is well known that different 
crystalline faces of a metal have distinct catalytic activ- 
ity. However, the present assumption is useful for a study 
which aims at investigating the interplay of many other 
different physico-chemical processes. 

III. ANALYTIC SOLUTION 

In order to solve the model analytically, we assume 
that the catalyst particles and the support segments be- 
tween them are sufficiently large (/ ^ 1, d ;2> 1), so that 
a continuous approximation is possible. Although I and 
d may not be very large in real systems, we will show 
that the continuous approximation works well even with 
I and d of order 10. 

The dimensionless reactant coverages in the catalyst 
particle and in the support are respectively defined as 
9c(x,t) and 9s{x,t). For simplicity, we use the same 
variable x for position in both regions, with the range 
—la/2 < X < la/2 in the catalyst and the range 
—da/2 < X < da/2 in the support. Diffusion, reaction 
and adsorption-desorption processes of Fig. 1 lead to 
equations for surface coverages; in the catalyst, we have 

^0c(a;, t) = D^e,{x, t)+F [1 - 9, (x, t)]-{k,+k'M^, 

(4) 

and in the support we have 

^^9six, t) = D-^9s (x, t) + F[l-9, {x, t)] - kl9, {x, t) . 

(5) 

Each of these equations is valid in the above defined 
ranges of x. 

Here we are interested in steady state solutions, where 
■^9c{x,t) — ■^9s{x,t) — and, consequently, 9c and 9s 
depend only on x. In the catalyst, this gives 

+ F{1^9c)- (k, + fc^)0c = . (6) 

dx 

In the support, an analogous equation (without the reac- 
tion term) is obtained. Eq. ^ can be easily solved and 
gives 

6'c(a;) = rc + Q;cCosh(a;/Ac), (7) 



where ac is a constant to be determined from boundary 
conditions and 

r 1 A ^ / 

" \ + h/F + kl/F ' ' \ \ + h/F + kl/F' 

(8) 

Analogously, in the support we obtain 

9s{x) ^rs + asCOsh{x/K), (9) 
where as is a constant and 

_ 1 / D/F 

l + ' ^^ = '^\ + kl/F- ^^^^ 

Note that, as expected, diffusion, reaction and desorp- 
tion rates appear in our results in the form of ratios to 
the external particle flux F. 

The calculation of unknown constants in equations 
such as d?]) and (O usually follows from the use of suitable 
boundary conditions. However, here diffusion leads to a 
net flux of reactants from the catalyst to the support or 
vice- versa, which depends on the competition of all other 
processes along both regions. Thus, those unknown con- 
stants will be determined by matching the gain and loss 
terms in each region due to all those processes, and the 
net flux at the interfaces will be obtained from them. 

In the catalyst, the contribution to the loss rate due 
to diffusion involves the probability of finding a reactant 
in the edge site of that region and of finding the neigh- 
boring support site empty (due to the excluded volume 
condition). Other contributions come from reaction and 
desorption along the particle. Thus, the loss in the cov- 
erage of the catalytic region per unit time is 

(Aec),_ = — (l-r)0t + _ / {k, + kl)9c{x')dx' 

« a J -la/2 

on 

= {i-9*)e^ + i{h + ki)9c. (11) 

where 9* is the coverage of the edge site of the support 
region (neighbor of the catalyst) 

9* ^9s{x^'^^^n + as cosh (^^) , (12) 

6^ is the coverage of the edge site of the catalytic region 
(neighbor of the support) 

et^^?c(^:r=|^ =rc + acCOsh(^^), (13) 

and 9c is the average coverage of the catalytic region 

1 pla/2 

Oc = ^ 9c{x')dx'. (14) 

l-O- J -la/2 

The gain rate, which accounts for flux from the support 
to the catalyst at the edge sites and for the external flux, 
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" " J -lall 



a 
2D 



f.la/2 
2/2 

'{l-0^)+lF{l-0c). 



(15) 



Analogously, loss and gain terms can be determined for 
the support region, so that solutions for ac and as are 



etc 



and 



Ac (r3-rc)tanh(^) 



A, cosh ( ^ ) tanh ( )+ [A^+a tanh ( ^ )] sinh ( ^ ) 



,(16) 



As(rc - rs)tanh ( ^ 



^= cosh ( ^ ) tanh ( ^ ) + [Ao +a tanh {^)] sinh ( ^ ) 

(17) 

The turnover frequency, which is the number of reac- 
tions per unit site and unit time, is given by 



TOF = ek^e^ 



(18) 



We simulated the discrete model, as defined in Sec. 2 
(see also Fig. 1), in order to check the accuracy of the 
analytic solution when I and d are not very large. Typical 
simulations consisted in 100 realizations of the process in 
a lattice with L = 2^^ sites and catalyst coverage near 
15%. At each step, an attempt to deposit a new reac- 
tant at a randomly chosen site is done. Subsequently, 
the numbers of attempts to move reactants, desorb them 
and perform reactions are chosen proportional to the re- 
spective rates (Z?, k^, kr). Simulations begin with 
an empty lattice and proceed up to a long time after a 
steady state has been reached, where the coverages of 
both regions are constant. 

We observe that even for / ^ 10 the continuous ap- 
proximation is good. This is illustrated in Fig. 2a and 
2b, where we compare the analytic and numerical results 
for the coverage distribution in the catalyst and in the 
support, respectively, using I = 18 and d = 100. Slight 
deviations are only found near the boundaries of those 
regions, but are always smaller than 5% for small I. The 
accuracy in the average coverages and in the TOF is usu- 
ally higher. 

Here the model was solved in the steady state, but it 
certainly can be extended to other situations, for instance 
when a finite amount of the rcactant flows to the catalyst 
surface and the TOF changes in time. Other possibili- 
ties are the explicit incorporation of sticking probabili- 
ties and the assumption of different diffusion coefficients 
in the particles and in the support. The one-dimensional 
structure facilitates the solution, while preserving essen- 
tial ingredients such as the spatial heterogeneity. 



IV. RESULTS 

In order to understand the role of the different rates 
from a scaling approach, we consider the diffusion lengths 
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FIG. 2: Reactant coverage as a function of the position x 
along (a) the catalyst and (b) the support, for 2D/ {a-^F^ = 
420, kl^/F = 0.1, kr/F = 8, d = 100 and I = 18. Circles are 
simulation results and the curves are analytical results. 



in the catalyst particle and in the support. If the corre- 
sponding region is large enough, the diffusion length mea- 
sures the typical distance a reactant moves on it before 
reacting or desorbing. On the other hand, if the length 
of that region is smaller than the diffusion length, then it 
is expected that the reactant reaches its border and can 
flow to a neighboring domain. 

In the catalyst, the average lifetime of an adsorbed 
species before reacting or desorbing is Tc ^ 1/ (fcr + fcj;)- 
During this time, it executes random walks with diffusion 
coefficient D, thus the average distance it spans is of 

1/2 

order {Dtc) . This is the so-called diffusion length. 



D 



kf, 



(19) 



Analogously, a reactant in the support has a typical life- 
time 1/A;^, thus the corresponding diffusion length is 



(20) 



These expressions can be easily generalized to the case 
of different diffusion coefficients in the particles and in 
the support, which means that interpretations based on 
these quantities have a broader applicability. 

We will consider cases where desorption in the catalytic 
region is very low, i. e. k^ <IC F^k^^krT^D/a^. The 
values of all rates presented below are given relatively to 
the incident flux rate F, thus they are all dimensionless 
(setting F = 1 and the value of the parameter a, we 
would get the other ones in SI units). The TOF is also 
expressed relatively to F, thus it is limited to a maximum 
TOF/F = 1. 

In order to facilitate the presentation of the results, 
hereafter we will use the terms spillover and back spillover 
to denote a net flux of reactants by diffusion from the 
catalyst to the support and vice-versa, respectively. As 
noted above, the broader use of these terms follows a 
trend of some recent works 0, [2l|, [H, [H, [H, ill . 
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A. Effects of reactant mobility and catalyst 
geometry 

First we distinguish the conditions where either 
spillover or back spillover is dominant as the reactant 
mobility in the surface increases. In Figs. 3a and 3b we 
show the normalized TOF as a function of 2D /a? for sev- 
eral reaction rates, respectively with fractions of support 
covered e = 0.05 and e = 0.15, and the same particle size 
I = 75. 
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FIG. 3: Normalized turnover frequency as a function of diffu- 
sion coefficient for several reaction rates, with fc^/F = 10"'^, 
k^^jF — 10 and I — 75. Fractions of support covered are (a) 
e = 0.05 and (b) e = 0.15. 



Even without reactant diffusion {D = 0), the TOF sig- 
nificantly varies with because a large area of the cata- 
lyst may be poisoned for low reaction rates; for instance, 
in Fig. 3a we have normalized TOF = 0.0278 and 0.0546 
for kr/F = 0.1 and kr/F = 5, respectively. Indeed, ex- 
cluded volume effects limit the adsorption process, which 
is known to be a relevant effect even in mean-field models 

. The increase of diffusion coefficient improves the cat- 
alytic process for high reaction rates, since back spillover 
effects are dominant as D increases. However, with low 
reaction rates, the opposite effect is observed: reactants 
more easily leave the particles as D increases, going to 
the support where they rapidly desorb. In Fig. 3a, there 
is no change in the normalized TOF as D increases for 
kr/F = 0.525, which corresponds to Lg/Lc ~ 0.23. If 
other values of the rates and other values of I are cho- 
sen, the same feature is observed for a different value of 
kr, but with the same ratio Lg/Lc- On the other hand, 
in Fig. 3b (e = 0.15), that feature is observed when 
Ls/Lc « 0.42. Thus, the critical ratio Ls/Lc which sepa- 
rates regimes of rapid and slow reactions depends only on 
the fraction of support covered e; above (below) the crit- 
ical ratio, back spillover (direct spillover) is dominant. 
This result can be interpreted as follows: if the rela- 
tive increase of diffusion length in the support is larger 
(smaller) than that in the catalyst, then more reactants 
flow towards the particles (support) and the TOF in- 
creases (decreases). 

Now we consider the effects of particle size. We con- 
sider changes in I with fixed e, in order to simulate cases 



where a fixed amount of catalytic material is deposited 
on the support, but islands of different sizes are formed. 
Under these conditions, it is important to stress that the 
size of the gaps between particles (d) also increases when 
I increases (Eq. [T]). 

In Figs. 4a and 4b we show the normalized TOF as a 
function of the particle size I for two different fractions 
of the support covered and various reaction rates. Again 
we observe regimes of high and low reaction rates, corre- 
sponding to high and low ratios Ls/Lc- For large reaction 
rates, decreasing I is favorable, since d also decreases (e 
is fixed) and facilitates the back spillover. On the other 
hand, for low reaction rates, the conversion is improved 
by increasing the particle size, since the reactants spend 
longer times in larger catalytic regions, which compen- 
sates the increased desorption in the support. In other 
words, the loss due to decreasing back spillover is com- 
pensated by a gain in decreasing spillover. The same 
ratios L^/Lc, which depend on the fraction of support 
covered e, separate the two regimes. 
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FIG. 4: Normalized turnover frequency as a function of cata- 
lyst particle size for several reaction rates, using k'^/F = 10~^, 
k^F = 1 and 2D/ (a^-P) = 500. Fractions of support cov- 
ered are (a) e = 0.15 and (b) e = 0.3. 

These nontrivial effects of the catalyst geometry may 
be helpful for catalyst design if one is interested in taking 
advantage of back spillover, preferrably in the case of 
rapid reactions, or in reducing the effect of spillover in 
the case of slow reactions. At this point, it is important 
to notice that Figs. 3 and 4 show that the change in one 
model parameter may increase the TOF by a factor near 
3, which is a remarkable change in the efficiency of the 
catalytic process. 

As far as we know, this is the first work that analyzes 
the conditions for changing the direction of the net reac- 
tant flux in the particle-support interface. The interpre- 
tation of results based on diffusion lengths was formerly 
proposed in Refs. [13, HI] for models of CO oxidation 
in oxide supports. However, both studies focused on the 
regime where back spillover is dominant, for instance by 
assuming infinite diffusion lengths of some species ad- 
sorbed in the particles. That regime was also considered 
in Ref. psj with a mean-field approach that accounts 
for the different neighborhood of the catalytic sites (in 
an approximate form). In this case, the approximation 
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is successful because the diffusion lengths are small. 



B. Temperature effects 

The diffusion coefficient and the reaction and desorp- 
tion rates are expected to have Arrhcnius forms as fol- 
lows: 



and 

Here, Vi is a frequency and Ei is an activation energy 
[i — D,d,r). Assuming that the activation energy for 
desorption in the catalyst particles is much larger than 
the other activation energies, we use k^F = 10"^ which 
is negligible compared to the other rates in the relevant 
temperature ranges. For the other activation energies, 
we always assume that Ed < Ed [ioj . 

Other reasonable assumptions on the amplitudes of 
the Arrhcnius rates facilitate the analysis of the effects 
of different ranges of energy barriers. We will consider 
vd = Vr = 2 X 10^^ ^ a = 5 A, and work with a range 
of ratios Vd/vr from 10 to 1000. These relations are rea- 
sonable for Langmuir-Hinshelwood (LH) reactions and 
for CO adsorption in oxides 0, [H, |4^, respectively, but 
we emphasize that these systems are only rough guides 
to choose working parameters, and may not be viewed as 
prospective applications - indeed, the unimolecular reac- 
tion of our model is representative of ER mechanism. In 
the following, we also consider fraction of support cov- 
ered e = 0.15 and F = 1 (thus the calculated rates 
are again ratios to F). 

Qualitatively, we expect that each microscopic pro- 
cess will be significantly activated when its rate exceeds 
the external flux F. However, excluded volume effects 
lead to surface poisoning when reactions are not frequent 
(low temperatures) , so that other processes can affect the 
turnover frequency only after reactions are activated. 



1. The cases Er > Ed and E,- < Ed 

First we consider the simplest case where reactions are 
very difficult compared to the other activated processes, 
i. e. Er > Ed- An example of slow reaction is the 
hydrogenation of CO on Pt , when compared to the 
spillover to the Ti02 support and the formation of CH^^O 
there. 



Fig. 5a shows the typical evolution of the TOE with 
temperature for different particle sizes. When the reac- 
tions become more frequent (fc^ ^ F, i. e. T ^ 270 K), 
desorption is already activated and diffusion is very fast. 
Thus spillover and subsequent desorption of reactants 
does not allow the increase of the turnover frequency. 
Instead, the TOE begins to increase only when kr '~ kd 
{T ^ 550 K). At higher temperatures, the negative con- 
tribution of spillover is more important when the par- 
ticles are small, so that a slow increase of the TOE is 
observed. On the other hand, this negative contribution 
is reduced for large particles, and the maximum TOE 
is rapidly attained by increasing the temperature. The 
value TOY / F w 0.15 corresponds to the fraction of the 
surface covered by the catalyst, which is expected at high 
temperatures because only species adsorbed on the cat- 
alyst react. 
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FIG. 5: Normalized turnover frequency as a function of tem- 
perature for several particle sizes, with: (a) Ed ~ 5 kcal/mol, 
Er = 15 kcal/mol and Ed = 10 kcal/mol; (b) Ed ~ 
5 kcal/mol, Er = 3 kcal/mol and Ed = 10 kcal/mol\ (c) 
Ed ~ 6 kcal/mol, Er = 3 kcal/mol and Ed = 35 kcal/mol. 
In all cases, Vd/^r = 100. 

Next we consider the case where reactions are easily 
excited, i. e. Er < Ejj. An example of rapid reaction 
with ER mechanism is CO oxidation on Pd/Ce02 (2l| . 
where reaction rates are nearly 100 times larger than back 
spillover rates of oxygen. 

Figs. 5b and 5c illustrate the case Er < Ed with small 
and large desorption energies, respectively. 

In Fig. 5b, the TOE increases towards a plateau in 
TOV/F « 0.15 at T ~ 50 - 60 K), where K ~ F. At 
T ^ 100 K, when diffusion is activated, back spillover 
leads to a second jump in the TOE, which is highly de- 
pendent on the particle size. For smaller sizes (small 
gaps between particles), the migration of reactants from 
the support to the catalyst is easy even for low D, thus 
large conversion rates are rapidly obtained. The tem- 
perature of maximal TOE is attained when desorption 
begins to play a significant role (fc^ ~ F), independently 
of particle size. In the case of large gaps between the par- 
ticles, this temperature is still low for back spillover to 
be efficient, thus only a small peak appears in the TOE 
plot. For further temperature increase, /-dependent re- 
sults are again obtained. The diffusion length in the sup- 
port is Ls = ^J^^exp[{Ed- Ed) /2kBT], which de- 
creases with increasing temperature, and the beneficial 
effect of back spillover ceases when ^ I; this condition 
is satistified at higher temperatures for smaller I, which 
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explains the slower decay of the TOF in this case. For 
these reasons, the peak in the TOF is high and broad for 
small /, low and narrow for large I. 

In Fig. 5c, the main features of Fig. 5b are present. 
However, since Ed is much larger than the other ac- 
tivation energies, the maximal effect of back spillover 
(TOF/F « 1) is observed in a wider temperature range 
and for all particle sizes shown there. Thus, Fig. 5b il- 
lustrates typical conditions in which particle size effects 
are clearer: desorption in the support has higher activa- 
tion energy than diffusion and reactions, but activation 
of a process begins while the other processes are not fully 
activated, so that the diffusion length Ls cannot attain 
large values before desorption is activated. 

2. The case of intermediate reaction energies 

Now we consider cases of intermediate activation en- 
ergies for reaction, i. e. Eu < Er < Ed- CO oxidation 
provides several examples with such relation between ac- 
tivation energies; however, it is important to stress the 
difference in the usual reaction mechanism (LH instead of 
ER). An example is CO oxidation on Pt/Ce02 [l9|: the 
energy of diffusion of O on the support is 18kcal/mol, 
while reaction energy is 27kcal/mol and desorption in 
the support is 60kcal/mol (diffusion on the catalyst is 
assumed to be fast for the application of a mean-field 
model in Ref. [29^). 

In Fig. 6a we show the normalized TOF as a function 
of temperature for our model with three different par- 
ticle sizes, using a set of activation energies previously 
suggested for CO oxidation in Pt{lll) and I'd — lOOi^r 
Fig. 6b shows the evolution of the average coverage 
(particle plus support). 




T/K T/K 

FIG. 6: (a) Normalized turnover frequency as a function of 
temperature for several particle sizes, with Ed ~ 6 kcal/raol, 
Er = 24 kcal/mol Ed = 35 kcal/mol, and Vdjvr = 100. (b) 
Temperature dependence of the corresponding average reac- 
tant coverage. 

At low temperatures, all the surface is poisoned, and 
reactions are slow. Increasing the temperature to T 
400 K (kr ~ i^), the turnover frequency increases, thus 
vacancies are created in the particles and in the support, 
which facilitates the back spillover effects. Note that 



there is no significant effect of particle size when a large 
temperature range is scanned, since the diffusion lengths 
Ls and L(. are very large (diffusion is highly activated 
in much lower temperatures). While the TOF increases 
with temperature, the coverage rapidly decreases towards 
zero because diffusion and reaction are rapid compared 
to the external flux. Further temperature increase leads 
to a maximum in the turnover frequency, again when des- 
orption in the support is activated (fc| ~ F\ However, in 
the right side of the peaks of Fig. 6a, we observe a size de- 
pendence because back spillover ceases only when Lc, ~ 
At high temperatures, a size dependence is also notice- 
able: small particles lose more reactants by spillover than 
the large ones, and those reactants easily desorb in the 
support, which leads to smaller TOF. 

As explained above, the most remarkable effects of par- 
ticle size are observed in cases where the activation en- 
ergies are close to each other. This is also illustrated in 
Figs. 7a and 7b for Eu < E^ < Ed, but Ed close to 
Er, with different ratios between i>r and i^d- In this case, 
when reactions become more frequent and leave room for 
new incoming reactants {k^ ~ F), the surface mobility 
is still low (i. e. IDjc? is not much larger than F\ 
Thus, back spillover is significant only for small parti- 
cles (and small vacancies between them). However, if 
the spacing between particles is large, then slow diffu- 
sion is not able to bring all reactants adsorbed in the 
support to the catalytic region. Figs. 7a and 7b also 
show that these effects appear in a wider temperature 
range when the amplitudes of the Boltzmann factors are 
closer to each other, since activation of desorption occurs 
in higher temperatures. 




T/K T/K 

FIG. 7: Normalized turnover frequency as a function of tem- 
perature for several particle sizes, with Eo — 5 kcal/mol, 
Er = 6 kcal/mol Ed ~ 10 kcal/mol: [a^Vd/i^r ~ 10; (b) 
Ud/l^r = 1000. 



V. CONCLUSION 

We proposed a model for a unimolecular reaction in 
a supported catalyst with tunable particle size and frac- 
tion of support covered. It represents important physico- 
chemical processes, such as diffusion, desorption and the 
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external flux of reactants. We analyzed the effects of 
spillover and back spillover on the turnover frequency 
under different conditions. 

First we considered the isolated effects of enhancing 
surface diffusion (coefficient D) and increasing catalyst 
particle size I. For rapid reactions, the increase of any of 
those quantities favors back spillover and, consequently, 
increases the turnover frequency. This regime is sepa- 
rated from that dominated by spillover by a critical value 
of diffusion length in the catalyst and in the support, 
which depends only on the fraction of the support sur- 
face covered by the catalyst. In the spillover-dominated 
regime (slow reactions) , increasing D or I slows down the 
conversion of the reactants. 

Subsequently, we considered temperature effects by 
assuming Arrhenius dependence of all physico-chemical 
rates and reasonable values for the amplitudes in those 
relations. With activation energy for desorption (Ed) in 
the support larger than that for diffusion (Eu), a peak 
in the turnover frequency as a function of temperature 
is observed for small and intermediate values of reaction 
activation energy E^., i. e. for cases where Er < Ed- Sig- 
nificant particle size dependence in the peaks is usually 
observed when the gaps between those energies are small, 
so that activation of one process occurs while the other 
ones are not fully activated yet, and the corresponding 
diffusion lengths rapidly vary with temperature. The 
right side of those peaks show size-dependence under 
more general conditions. For fixed amount of catalytic 
material deposited on the support, small particle sizes 
(with small distance between them) allow the turnover 
frequency to attain high peak values due to the benefi- 
cial effect of back spillover, while large particles provide 
low enhancements of catalytic activity. Finally, in the 
case of slow reactions [E^ > Ed), the TOF monotoni- 
cally increases with temperature, and large particle sizes 
are more efficient to avoid the negative effects of direct 
spillover. 

The aim of the present work is not the quantitative de- 
scription of a particular catalytic process, but to discuss 
the interplay of various physico-chemical mechanisms 
and conditions in systems where spillover is present. 
However, it may be useful in the interpretation of some 
experimental results and motivate the proposal of ex- 
tended models for their quantitative description. 

Here we mention two recent works as examples where 
the scaling ideas developed above may be relevant. How- 
ever, we emphasize the fact that our model cannot be 
directly applied to these systems, thus our results only 
suggest general guidelines to understand their qualita- 
tive behavior. In the first example. Piccolo and Henry 
studied the oxidation of CO by NO on Pd/MgO and 
observed a peak in the TOF as a function of the tem- 



perature. They found a remarkable increase in the peak 
height as the particle size was decreased (even being ac- 
companied by a decrease in the Pd coverage), but tiny 
shifts in the temperature of the maximum. These results 
resemble those in Figs. 6a, 7a and 7b. While the size- 
independent curve shape suggests weak size-dependence 
of activation energies, the peak height increase suggests 
that diffusion lengths of reactants were not very large 
when reactions were activated, so that back spillover 
(which increases TOF for any type of reaction) was facili- 
tated with small particles. Thus, differences in activation 
energies of diffusion, reaction and desorption are proba- 
bly small, with Ed being larger than the other ones. The 
second example is a study of CO oxidation on Au on 
active supports (i. e. those which adsorb and supply 
oxygen to the reaction), which shows a different trend 
[47| : under certain conditions, the TOF does not depend 
on catalyst particle size. Although that reaction involves 
two species and is probably of LH type, the experimen- 
tal result suggests that diffusion was highly activated at 
the working temperature. In this situation, similarly to 
our model, diffusion lengths are always very large and the 
contribution of back spillover does not depend on particle 
size. 

Recalling the results of models for ethene and acetylene 
hydrogenation which in corp orate spillover effects is also 
interesting at this point [35[3l|. In both systems, the hy- 
drocarbon blocks catalytic sites at high pressures, which 
leads to a discontinuity of the TOF when the external 
reactant flux towards a catalytic site is of the same order 
of the microscopic reaction rate. In ethene hydrogena- 
tion, this effect is shown to be more intense with small 
particles [SO], leading to finite-size effects. However, dif- 
fusion is weakly activated at the working temperature 
when compared to the other processes (the amplitudes 
of Boltzmann factors associated to different processes are 
very different in that case), thus there is no significant 
effect of diffusion in the TOF. This reinforces the conclu- 
sion that the scaling approach proposed here can be ex- 
tended to interpret more complex reaction models. This 
is very important in cases where analytic solutions are 
not feasible because it may help choosing the conditions 
to perform simulation work. 
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